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Abstract 



Motivated by the need for an absolute polarimeter to determine the beam polar- 
ization for the forthcoming RHIC spin program, we study the spin dependence of the 
proton-proton elastic scattering amplitudes at high energy and small momentum trans- 
fer. In particular, we examine experimental evidence for the existence of an asymptotic 
part of the helicity-flip amplitude ^5 which is not negligible relative to the largely imag- 
inary average non-flip amplitude cf> + = \{(f>i + 4>z)- We discuss theoretical estimates 

of rs = m (p5 jyf—i Im cp + based upon several approaches: extrapolation of low and 
medium energy Regge phenomenological results to high energies, models based on a 
hybrid of perturbative QCD and non-relativistic quark models, and models based on 
eikonalization techniques. We also apply the rigorous, model-independent methods of 
analyticity and unitarity. We find the preponderence of evidence at currently avail- 
able energy indicates that rs is small, probably less than 10%. The best available 
experimental limit comes from Fermilab E704: combined with rather weak theoretical 
assumptions those data indicate that \r^\ < 15%. These bounds are important because 
rigorous methods allow much larger values. Furthermore, in contradiction to a widely- 
held prejudice that r§ decreases with energy, general principles allow it to grow as fast 
as In s asymptotically, and some of the models we consider show an even faster growth 
in the RHIC range. One needs a more precise measurement of rs or to bound it to be 
smaller than 5% in order to use the classical Coulomb-nuclear interference technique 
for RHIC polarimetry. Our results show how important the measurements of spin 
dependence at RHIC will be to our understanding of proton structure and scattering 
dynamics. As part of this study, we demonstrate the surprising result that proton- 
proton elastic scattering is self-analysing, in the sense that all the helicity amplitudes 
can, in principle, be determined experimentally at small momentum transfer without 
a knowledge of the magnitude of the beam and target polarization. 
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1 Introduction 



The need to understand the spin dependence of scattering amplitudes at high energy and 
small momentum transfer is important for two distinct reasons. Firstly it is a great challenge 
to strong interaction theory, since it involves the application of QCD in a kinematical region 
where non-perturbative effects are important. QCD has had great success in the perturbative 
region, but experiments at HERA at very small x are already raising questions for which 
the standard perturbative approach may be inadequate Q; and future experiments at RHIC 
and LHC will produce a vast amount of data outside the perturbative region. It is hard to 
imagine a global solution to a non-perturbative QCD effect such as small-t spin dependence, 
but it is becoming more and more urgent to try to make some progress in this direction. 

Secondly the extremely important RHIC spin program [fj, [J, which will test many 
elements of QCD at a new level of accuracy and detail, relies heavily upon an accurate knowl- 
edge of the beam polarization. For the purpose of measuring the beam polarization P, the 
Coulomb-Nuclear Interference (CNI) polarimeter is very attractive: it has a reasonably large 
analyzing power (about 4%) in a region of momentum transfer (|i| « 0.002 — 0.003 GeV 2 ) 
where the rate is extremely high. This method depends on the dominance of the interfer- 
ence of the one-photon exchange helicity-flip amplitude (by an abuse of the term, normally 
called a Coulomb amplitude, more properly the magnetic amplitude) with the non-flip strong 
hadronic amplitude, which is determined by the total cross section. The accuracy of the 
method is limited by our uncertain knowledge of the hadronic helicity-flip amplitude; its 
interference with the non-flip one-photon exchange amplitude has the same shape in this t- 
region 0, ^, || [5], and so must be known, or limited in size, in order to achieve the required 
accuracy. The requirements of RHIC polarimetry (AP/P < 0.05) §J put very stringent de- 
mands on our knowledge of the helicity-flip amplitude. This problem was the impetus that 
drew our attention to the long-standing question of the size of the proton-proton helicity-flip 
amplitudes. 

This is not intended to be a paper on polarimetry, though we will inevitably make 
further comments on the subject as appropriate; indeed, the demands of RHIC just cited set 
a standard for our investigation. The aim of the paper is to provide a reliable assessment 
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of what is known about the helicity-flip amplitudes and what is expected for them at high 
energies on the basis of various approximate or rigorous theoretical calculations. 

Another well-known practical issue arising from our lack of knowledge of spin de- 
pendence is in the determination of the total cross section via the use of unitarity and the 
extrapolation of the differential cross section 0, [10], [11]. In particular, this may lead to an 
overestimate of the total cross section by an amount proportional to the ratio of the sum of 
the squares of the helicity-flip amplitudes to the square of the non-flip amplitude at t — 0. 
To put this statement more correctly and more precisely, in well-known notation which will 



be fully defined in Section 2, it will be overestimated by the factor [T2| 

Jt + P 1 (i) 

where 

B 2 = 1 /AaL Wl+P 2 -) 1 /A^A 2 (1 + pf) 

P 4UJ (1 + p 2 ) ^2 U to J (1 + p 2 )' 1 ' 

Martin 0] has emphasized that, because this is a ratio of squares, a quite good comparison 
between cross-sections obtained by this technique and more direct measurements of a tot 
leaves room for substantial spin dependence. 

Both of these experimental issues along with the theoretical studies using unitarity 
and dispersion relations emphasize the importance of understanding spin dependence at very 
small \t\. In addition, the very powerful tool of interference between Coulomb and strong 
amplitudes for extracting small parameters (like the p parameter for unpolarized elastic 
scattering) is effective in this region. 

Of course the interest of this physics has been understood for a very long time. The 
earliest studies relevant to our work date from the sixties. Associated in large part with the 
polarized proton programs at Argonne, CERN and Serphukov, there was a very large amount 
of phenomenological work in the seventies, and there were at the same time a number of new, 
fundamental ideas introduced. In the eighties and later, QCD has led to new techniques for 
modeling the spin dependence of high energy scattering, and the experimental program at 
Fermilab has made important contributions in this field. Specific citations will be given at 
the appropriate place in the following sections. With the coming of RHIC, the experimental 



motivation is very strong to revisit past studies and to attempt to make some advances on 
them. That is our purpose here. 

Section 2 will lay the groundwork for subsequent discussion by defining the basic 
amplitudes and expressing the various measurable polarization dependent quantities in terms 
of them. The general forms near t = will be discussed using Regge concepts, especially 
charge conjugation C and signature (— 1) J of the exchanged system, and the implications 
for the asymptotic phase of the various amplitudes. The terms "pomeron" and "froissaron" 
will be defined for our purposes, and several general results will be reviewed. 

In Section 3 our best knowledge regarding helicity-flip amplitudes will be given. This 
includes low and moderate energy Regge and amplitude analysis for pp and 7rp scattering, 
the energy dependence of P = An at small t and the most pertinent piece of experimental 
information: the measurement by E704 at Fermilab of in the CNI region. 

Section 4 applies the rigorous methods used to derive the Froissart-Martin bound 
to limit the energy dependence of the single helicity-flip amplitude relevant for CNI, and 
interprets this in terms of the impact parameter representation. 

Section 5 contains a description and evaluation of several models which give predic- 
tions for spin dependence at high energy. These will mainly address the single helicity-flip 
amplitude relevant for the CNI polarimetry. 

Section 6 reviews the issues of Coulomb enhancement and shows how, in principle, 
all the scattering amplitudes in pp scattering may be determined experimentally without 
knowledge of the beam polarization P. This method is contingent on being able to make 
measurements of very likely tiny asymmetries and it may turn out not to be practical. 
Should such determination prove to be practical, elastic pp scattering could be used as a 
self-calibrating polarimeter. 

Lastly Section 7 gives our conclusions. 

Before moving on to the body of the paper, we would like to say that this work 
originated at a workshop sponsored by the RIKEN BNL Research Center during the summer 



of 1997 [fL4 |. During the workshop we, along with several other people, discussed and analyzed 
various other methods of polarimetry. Some methods are very clean theoretically and have 
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good analyzing power; in particular, polarized hydrogen jet targets provide a self-calibrating 
method |15), while elastic ep scattering is calculable and has a very large analyzing power 



with longitudinal polarized electrons and transverse or longitudinally polarized protons fl6 |. 
One can also calibrate an unpolarized hydrogen target with a second low energy scattering 
off Carbon; this requires working at larger \t\ where the rate is much lower, but values of 
t for which the analyzing power is large are sure to exist, in particular in the dip region 



[17]. Nuclear targets, either in colliding beam or fixed target modes, might be useful for 
elastic scattering in the same way, using structure at larger t; their use in the CNI region is 



subject to the same uncertainties as for pp [|17], [18], [19|| . Finally, because the purely empirical 
asymmetry observed in inclusive ir production is very large and the rate is high, it may be the 
most practical initial polarimeter f20|| ; it nearly meets the required precision standard but 
one needs data to calibrate this polarimeter using the same target and at the same energy 
(in the fixed target mode) as will be used in RHIC. The choice of method obviously involves 
several different kinds of factors some of which, such as technical and cost, are beyond the 
scope of this paper. 



2 Fundamentals and dynamical mechanisms 

It has long been understood that the measurement of helicity amplitudes at high energy 
could be a powerful tool for determining the dynamical mechanisms for scattering in the 
asymptotic region |2T], |22], |23|; this is especially true for nucleon-nucleon scattering because 
its very rich spin structure allows for a greater variety of quantum numbers to be exchanged 



24]. Five independent helicity amplitudes are required to describe proton-proton elastic 



scattering [|, 



i(M) = 


(++ \M\ ++), 


2 (M) = 


(++ \M\ --), 


s(M) = 


(+- |M| + -), 


4(M) = 


(+- \M\ -+), 


b(M) = 


(++ |M| + -). 



(3) 
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Here we use the normalization of H]. Since we are interested only in very high energy s/s, 
such as will be available at RHIC, and very small momentum transfer \t\ < 0.05 GeV 2 , we 
will generally neglect m with respect to s and neglect t with respect to m to simplify the 



presentation of the formulas which follow. For example, k 2 = J s(s — 4m 2 )/ '4 will be replaced 
by s/4. Then 

a tot = — lm(0i(s,t) + fo(s,t))\ t=Q (4) 



and 



% = 2 ^m 2 + \h* + \h 2 + n 2 + 4 w 2 }. (5) 



We will also have occasion to discuss (z) scattering of unlike- fermions, requiring a sixth 
amplitude 06, a single helicity-flip amplitude which degenerates to —05 for identical particles 
(of course, pp elastic scattering requires only 5 amplitudes), and (m) scattering of a proton 
on a spin-zero particle, like a pion or a spinless nucleus, requiring only two amplitudes, a 
non-flip and a flip amplitude. 

We will consider only initial state polarization measurements. There are certainly 
interesting things that can be said about final state polarizations, but the first generation 
spin program at RHIC will not measure these and so we will not discuss them here. Using 
only initial state polarization, with one or both beams polarized, one can measure seven spin 
dependent asymmetries. We follow the notation of ||. There are slight variations in the 
definitions used in the literature, having to do with the orientation of axes. 

d(J 4:71 

AN ~dt = -^ Im {05(01+02 + 03-04)}, 
(1(1 A-7T 

A NN — = ^{2|0 5 | 2 + Re(0t0 2 -0^ 4 )}, 



Ass^- = J Re{0i02 + 0301}, 



d(J A.7T 
ASL ltt = ^ Re ^(^l + 02-03 + 04)}, 



All^ = ^{|0ir + |0 2 | 2 -|03| 2 -|04| 2 }. ( 6 ) 
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It will be convenient to introduce some shorthand: 

+ = i(0 1 + 3 ), 0_ = 1(^-03), (7) 

and 

_ Re 2 _ Re0„ 

^ 2 Im 02 ' ^ Im 0_ 

There are also two cross section differences corresponding to longitudinal and transverse 
polarization: 

Im0_O,O) lAdJs) 



Im0 + (s,O) 2 cr tot (s) 
Im0 2 (s,O) Acr T (s) 



At l =(x—^, (9) 



Aa T = a n -a n . (10) 



Im0 + (s, 0) (XtotO) 

When the proton scatters elastically off a distinct spin 1/2 particle, there are two more 
measurable asymmetries: A' N and Als, in obvious notation; these degenerate into A^ and 
Asl respectively when the two particles are identical. For scattering off a spin zero particle, 
there is only one asymmetry which corresponds to Ajq. 

At these small values of t, the interference of the strong amplitudes with the single 
photon exchange amplitudes will be important; this interference is central to this paper. To 
lowest order in a, the fine structure constant, one replaces 

l ->0 l + 0f n exp(z5) (11) 

with hadronic and electromagnetic elements. The Coulomb phase 5 is approximately inde- 



pendent of helicity [pi 26 



6 = aln q >(B + 8/A*) - a7 (12) 

where B, often called "the slope", is the logarithmic derivative of the differential cross section 
at t — 0, a number about 13 GeV -2 and increasing through the RHIC region, q 2 = —t, Euler's 
constant 7 = 0.5772 . . . and A 2 = 0.71 GeV 2 reproduces the small momentum transfer 
dependence of the proton form factors assumed to satisfy 

G E (q 2 ) = G M (q 2 )/v P = (1 + q 2 /A 2 )- 2 - (13) 



For pp scattering at high s and small t, the electromagnetic amplitudes are approximately 



rv c 

^ m = 0^ m = _F X 2 , 



2 

/em /em aSK 



Ft 



'2 ^ ~ 4m 2^2, 

asK 



FiF 2 , (14) 



where fi p — k+ 1 is the proton's magnetic moment, and m its mass. For the full expressions 
see, e.g., ||. The proton electromagnetic form factors F\(q 2 ) and i^*? 2 ) are related to Gg 



and Gm 0' sec ti° n 12.2] by 



j-, _ Ge — Gm t/ Am 2 _ Gm — Ge n _s 

Fl " l-t/4m* ' KF2 ~ l-t/Am 2 ' (15) 

The relations between 0! and 03 and between 2 and 04, Eq. (fHj), are special consequences 
of the quantum numbers of the exchanged photon; they are not generally true for the full 
amplitudes. Relations of this type will be dealt with shortly. 

Each hadronic amplitude 0-, can, in principle, be broken up into two parts 

0^0f + 0f (16) 

where <pf is controlled by Regge pole type dynamics and, in our normalization, decreases with 
energy roughly like s -1 / 2 with respect to the asymptotic part <pf s . Although the first term 
is essential to understanding the data in the low-to-moderate energy region which overlaps 
the RHIC range, we will focus here solely on the second term. 

Consider first the dominant non-flip forward amplitude + ; this must have an asymp- 
totic piece whose imaginary part grows with energy as a consequence of its connection Eq. (^) 
to the nucleon-nucleon total cross section. There are two widely used forms for 0^ s to de- 
scribe the high energy behavior of cr to t(pp), which is flat up to s/s ~ 20 GeV, with a value 
of 38 mb and then grows to 43 mb at y/s = 63 GeV increasing further to about 62 mb at 
the CERN SppS collider (y/s = 546 GeV). In the first, one fits the data with terms of the 
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form s ln p s, p < 2 |28], |23|. This form is suggested by Regge theory and the Froissart-Martin 
bound [0 

|0 + |<csln 2 s as s — > oo. (17) 

In this approach Im0+ s receives contributions from the simple pomeron pole IP, with inter- 
cept ap(0) = 1, together with a contribution growing at the maximum allowed rate sln 2 s 
(sometimes referred to as a froissaron [pgfl ) 

Im0+ s (s) = a P s + ap sin 2 s. (18) 



In the second, one introduces an "effective" pole, the Landshoff-Donnachie pomeron |3T| , 
with «p = 1 + Af>, where typically ~ 0.08. The ensuing behavior 

lm^ s ocs 1+Ap (19) 

gives an excellent description of the behavior of cr tot {pp) and cr tot (pp) and many other reac- 
tions. This form is also suggested by perturbative QCD calculations [j32|| , but with a larger 
value of Ap. However, ultimately, it violates Eq. fll7|) and so must be modified at higher 
values of s. This sort of behavior was obtained much earlier in QED-like theories [33j where 
consistency with Eq. fll7D was achieved through eikonalizing the form Eq. (|i"9f) . The unita- 
rization by multi-pomeron exchange of a "bare" pomeron which grows as s 1+Ap , Ap > 0, is 
obtained by eikonal methods in |}4, in those papers the relation of this result, via unitar- 



ity, to multiplicity distributions and inclusive inelastic cross sections is demonstrated. The 
resulting behavior is consistent with the Froissart-Martin bound, Eq. ([T7|) but the approach 
to the limiting asymptotic form is much more complex than is assumed in Eq. (|T%D. See the 
discussion later in Sections 4 and 5 and references cited there regarding the eikonalization 
method. 

There is also theoretical evidence, from a study of three-gluon exchange in QCD |36|, 



for a crossing-odd contribution to <p^ s which grows with energy slightly less rapidly than the 

pomeron exchange, and which would lead to a very slow decrease of the quantity 

(fltot(pp) — cr tot(pp))/(o"tot(pp) + o"tot(pp)) at asymptotic energies. However, phenomenological 



studies of this so-called odderon O contribution [29, [37| suggest that in the RHIC energy 
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range its contribution is very small compared to the crossing-even part of 0+ s . Roughly 

| ^pAs | odd 



< 2% (20) 



in the RHIC region and we shall therefore neglect the crossing-odd contribution to + in 
what follows. 

The key question for us is, do any of the non-dominant amplitudes 2 , <f>- and, 
especially, 5 have asymptotic behavior characteristic of the pomeron or froissaron? There 
is abundant evidence at low energy, some of which we will discuss in Section 3, that these 
amplitudes fall off with energy with respect to + as one would expect from lower lying 
Regge-exchange. It is not known, however, whether asymptotically they have a small but 
non-zero ratio to <f> + . To characterize these amplitudes we will define relative amplitudes in 
the following way: 



r 2 = R 2 + %h = 
r_ = R_ + il_ = 
r$ = Rb + ih = 



>2 



2 Im 0+ 



Im + 



m 0k 



t Im0 + ' 



r 4 =R 4 + i I 4 = (21) 

t lm0 + 

Notice the factor 2 in the definition of r 2 which is there to simplify many later formulas. 
The factors involving t which have been extracted reflect the fact that as t — * the strong 
amplitudes 0i, <p 2 and 03 go to a possibly non-zero constant while 04 oc t and 05 oc as 
a consequence of angular momentum conservation. The various r's will be assumed to be 
complex and to vary with energy but their variation with t over the small region we consider 
will usually be neglected. See, however, Section 6. 

The determination of the asymptotic spin dependence can be used to help identify the 
dynamical mechanisms at work at high energy. We can classify the dynamical mechanisms 
according to the the quantum numbers parity (P), charge conjugation (C) and signature (r) 
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of the t-channel exchange. An amplitude A T is called even or odd under crossing according 
as r = +1 or —1, since 

A T (e i7T s,t) = rA*(s,t). (22) 



For nucleon-nucleon scattering there are three classes of exchanges |23], |38[ and they 
contribute to the amplitudes as shown in Table 1. 



Class 1 

T = P = C 


Class 2 
T = -P = -C 


Class 3 
T = -P = C 


0+,05,02~04 


0- 


02 + 04 


P,0,p,u,f, a 2 


CLi 


■K,rj,b 



Table 1: Classification of pp amplitudes by exchange symmetries and the associated Regge 
poles 

If the asymptotically dominant contribution has definite quantum numbers, then unitarity 



requires that it has the quantum numbers of the vacuum |39| ; this is the defining property 
of the pomeron. Note that it is the quantum number C which determines the relative sign 
of the contribution of a given exchange to nucleon-antinucleon scattering i.e. 



(23) 



This implies that pomeron dominance and the absence of an odderon requires not only 
that the total cross sections for pp and pp be equal, but also their real parts, or p values. 
Because the pomeron has r = +1, the well-known argument relating the phase of a scat- 
tering amplitude to its energy dependence, see e.g. H40| , tells us that, if the asymptotic 



behavior of (a pp + a Pp ) goes as s a ln p s, then the amplitude for C — +1 exchange goes as 
s a ln p s exp(— tan/2)(l —tp7c/2\ns). Either of the two behaviors Eq. ([19]) or Eq. ([18]) imply 
that at the maximum RHIC energy range 



Re. 



P 



Im 



kAs 



0.12, 



(24) 



but the energy dependence over the entire range is somewhat different. (Of course, a detailed 
fit over the entire RHIC range will require the inclusion of lower lying Regge trajectories.) 
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It is not known whether the pomeron couples to 05 or to 02 — 04- The phenomeno- 
logical success at medium energies of "s-channel helicity conservation" |41| would suggest 
a small coupling, but this question is open to experimental study. If they do couple to the 
pomeron they will have exactly the same asymptotic phase as + . This may prove useful 
in investigating whether or not the dominant behavior becomes pure pomeron/froissaron as 
s — > oo, or if there can be substantial odderon contribution to these subdominant ampli- 
tudes. An odderon with nearly the same asymptotic behavior as the pomeron/froissaron 
will be approximately 7r/2 out of phase with it. As we have noted its coupling to + is 
quite weak, but nothing at all is known about its coupling to 02 — 04 or 5 and these phase 
relations may prove useful in probing for such couplings. This matter is of great interest and 
is discussed in a separate paper . 



The exchanged objects with the quantum numbers assignments in Table 1 could be 
pure Regge poles or cuts generated by the exchange of the Regge pole plus any number of 
pomerons. These cuts will have an asymptotic behavior which differs only by a power of 
Ins from the simple Regge pole and so must be considered along with it jKJ. In general, 
although the couplings of pure poles factorize, there is no reason for the cut couplings to 
do so. It is obvious that the charge conjugation parity of a cut is equal to the product of 
that of the poles that produce it. The corresponding situation with signature and parity is 
less obvious because of the relative orbital angular momentum the exchanged poles can have 
ZTJ. It has been shown, however, [4^, |46| that the signature of the cut r cut = T pole . This 



means that the important relation between C and r, that distinguishes Classes 1 and 3 from 
Class 2 in Table 1, is preserved for the cuts. The situation for parity is not as certain; Jones 
and Landshoff [47 have shown that the "wrong" parity cut, P cut = —P po i e is suppressed 
compared to the "right" parity cut, P cut = +P po i e - The strength of the suppression remains 
a quantitative question which is open to experimental and theoretical study. 

There are some very general things one can say about how the spin dependence can 
help distinguish pole from cut contributions; for an early example see If factorization 
should hold to a good approximation then one has 

2 M) = -|^4t and 0-M) = O. (25) 
<p+{s, tj 
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This obviously leads to a very simple spin dependence. In particular it implies, that as t — > 0, 
4>2 oct rather than the generally allowed behavior. 

Even if factorization is not valid, some of the same conclusions can be obtained just 
on the basis of quantum numbers. One particularly important example has to do with 02 
and 04. We have had little to say about 04 because angular momentum conservation forces it 
to vanish linearly as t — > 0. If either factorization holds or the dominant exchange has pure 
CP = 1 or CP = —1, then 2 must also vanish in the forward direction [32, The first 



condition we have just seen. The second can be confirmed by examining the table. There one 
sees that 2 + 04 and 02 — 04 couple to opposite values of CP. Therefore if only one value of 
CP is dominant asymptotically, 02 ~ =F04 as s — > oo and it, too, must vanish at t = 0. This 
makes the measurement of 02 near i = 0a very interesting probe of the dynamics; it may, at 
the same time have the unfortunate side effect of making some asymmetries unmeasurably 
small. 

Finally, notice from the table that neither the pomeron nor the odderon have the 
quantum numbers required to couple to 0_; it thus seems unavoidable that 

IQti 

Aa L = Im0_ (26) 

s 

should vanish like s _1//2 as s — > oo. This we have seen is also a consequence of factorization 
p2fl . If it does not, it indicates an asymptotically important exchange other than the pomeron 
or the odderon. Such an object has never been suggested to our knowledge, but there is no 
obvious reason that it should not exist. 

We see here some very simple statements that one can make which characterize the 
dynamics of high energy scattering by means of the spin variables. If the dynamics is well 
approximated by a pure pomeron pole the spin asymmetries will be quite small and require 
very sensitive experiments to measure. One should note that various suppressions as in 



pomeron vs. odderon or pole vs. cut |47] become gradually stronger (logarithmically with s 
or as a very small power of s); it will therefore be important to make these measurements 
over as wide an energy range as possible. RHIC presents a wonderful opportunity to do this. 
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3 Best experimental knowledge of 0_, 02 and 05 



As we have seen above, all the various spin observables are expressed in terms of the helicity- 
flip amplitudes. Clearly, to achieve a full amplitude analysis, one needs a substantial number 
of measurements, in the same kinematic region which is, unfortunately, far from the present 
experimental situation. Nevertheless, it is possible to extract from the available data some 
very useful information on the helicity-flip amplitudes which we will now try to review and 
summarize. 

Among the different spin observables we will consider, the transverse single-spin asym- 
metry An (or "analyzing power") has been extensively measured for pp elastic scattering, 
so it will play a central role in the following discussion. 



3.1 An in the CNI region 



The only experiment which has obtained relevant data in this kinematic region 
where —t is around 3 x 10~ 3 GeV 2 , is E704 at Fermilab ||8| at a lab momentum 



Pl = 200 GeV/c; the results are shown in Fig.l, along with two curves which will be ex- 
plained shortly. The errors are unfortunately too large to allow an unambiguous theoretical 



.06 



A 



N 0.04 



0.02 




0.01 



0.02 



0.03 



0.04 



0.05 



Figure 1: The data points are from Fermilab E704 pq/ . The solid curve is the best fit with 
the hadronic amplitude 05 constrained to be in phase with hadronic <p + ; the dotted curve 
is the best fit without this constraint. 
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interpretation, but let us now briefly recall what can one learn from it. From the formulae 
in Section 2, An is given by the expression (this is identical to the expression for the final 
state polarization parameter P) 



Im{(20 + + 2e^07 + 2 )* (0 5 + e^0 e 5 m )} (2?) 
N |0 + + e i5 e + m | 2 + l0-l 2 + |l02| 2 + 2|0 5 + e i5 0i m | 2 ' 1 ' 

for not too large values of —t, such that the amplitude 04 = (H — |0| — h) may be ignored 
because of the kinematical factor (— t) occuring in this double helicity-flip amplitude. 

In the one-photon exchange approximation 0^ m and 0| m are real and have well estab- 



lished expressions Eq. (14), so in order to make a theoretical prediction using Eq. (27), one 
needs to know the hadronic amplitudes + , 0_, 2 and 05. The imaginary part of the largest 
one, + , is related at t = to the total cross section <7 to t and the interference between 0| m 
and 0+ is most prominent when t = t c , where t c = — 9>ixa/a tQt . 

The explicit expression can be obtained by substituting the expressions from Section 
2 into Eq. (g7|): 

mAN 167T da R . r . . rT ^ , „ , T rT ^ ., t c 

— — — e~ Bt = [/t(l-<5p + Imr 2 -<mer 2 )-2(Imr 5 -<5Rer 5 )]- 
t c"tot dt t 

2(1 + Imr 2 ) Rers + 2(p + Rer 2 ) Imrs, 

2{p + 5) t i+ (l +p 2 )(l + /3 2 ), (28) 



where /3 is defined in Eq. (Q) . The asymmetry for the CNI region can thus be expressed 



as a quotient of a linear expression in t c /t in the numerator and a quadratic expression for 
t c /t in the denominator, neglecting terms of order t. 

The Coulomb phase 5 is small, about 0.02 in the CNI region, smaller at larger \t\. It 
has a slight effect on the position of the maximum in A^: 

tma: 



tr 



v / 3 + -(pImr 5 -Rer 5 )-(p + 5), (29) 



in the approximation where small quantities are kept to first order, but it enters the numer- 
ator multiplied by small amplitudes and so can be neglected for pp scattering. The height 
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of the peak is mainly sensitive to the unknown quantities lmr2 and Imrs, while the shape 
depends mainly on Rer 5 . For example, an Imr 5 value of ±0.1 modifies the maximum of 
by about 11%. 

There are two fits to the E704 data allowing a non-zero r^ shown in Fig.l M; the 
other rj's are set to zero. The solid curve is the best fit subject to the constraint that 05 is in 
phase with <f> + . The arguments in Section 2 show that if + and 05 have the same asymptotic 
behavior they will have the same phase; in that case the best fit is |r 5 | = 0.0 ± 0.16. Fitting 
without that constraint yields the dotted curve, which corresponds |r 5 | = 0.2 ± 0.3 with 
a relative phase angle to + of 0.15 ± 0.27 radians. Note the large uncertainties on these 
values. This is essentially the same as an earlier fit obtained in [[50] . As emphasized in |7], |§ , 



we see that a large value of Imr 5 generates a very large uncertainty on A max , which can be 
of the order of 30% or more. 

3.2 Energy dependence of the spin flip amplitudes from nucleon-nucleon scattering 

In the small t region we have some miscellaneous data on their magnitude and energy 
dependence. First, the transverse-spin total cross sections difference Acxt is related to Imr2 
for t — 0, according to Imr2 = — Aa^/2a to t- From the limited ZGS data [5T|, we find 



that Imr 2 decreases strongly in magnitude from —6% at pi = 2 GeV/c to —0.4% at pi = 
6 GeV/c. One can speculate whether for higher energy, it will remain negative and small 
or change sign and increase in magnitude. The charge exchange reaction np — > pn, can be 
also used to evaluate the modulus of 02 which dominates the cross section near the forward 
direction. The analysis of the data leads to the value \r 2 \ = 3.5% at pi, = 25 GeV/c 



and | 1 = 0.6% at pi = 270 GeV/c, further evidence for a strong energy fall off of the 1 = 1 
exchange amplitude. 

The longitudinal- spin total cross sections difference Ao"l is related to Imr_ for t = 0, 
according to Imr_ = Aa\ J /2atot- From the ZGS data [pj, we find that Imr_ decreases 
strongly in magnitude, from —10% at p^ = 2 GeV/c to —0.6% at pl = 12 GeV/c. At higher 
energy, E704 has measured A<tl pl for both pp and pp; their values imply that Imr_ has 



decreased below 10 3 for pp and to about 10 3 (with a 100 % error) for pp. These findings 
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Figure 2: An = P in percent as a function of lab momentum pi at t = —0.15 GeV 2 

are consistent with the belief that Im0_ vanishes as s — > oo. 

Away from the forward direction and the CNI region, the data indicate that An in pp 
elastic scattering is falling very fast with energy. This has sometimes led to the conclusion 
that the helicity-flip amplitude 05 would vanish as a power of s as s — > oo [[55]. In order to 
investigate this, we have taken a collection of data from various experiments which measure 
P = An at different energies, all for t = —0.15 GeV 2 (or interpolated from nearby values), 
the smallest \t\ for which there is sufficient data to do this p6fl . We have tried a fit suggested 
by Regge poles, namely P = a + b/ y/pZ + c/pl [0. This is shown in Fig. 2 and the relevant 
result is that a = 0.023 ± 0.012. It is not very well determined: it is consistent with pure 
CNI, which is approximately 0.01 at this value of t and pi ~ 300 GeV/c. At the same time 
it is consistent with a very large hadronic helicity-flip amplitude: the calculated value of 
An with Imr 5 = —0.6 and Rer 5 = —.015 (so that 05 is in phase with X ) approximates 
the fit very well for pi above 200 GeV/c . Because of the phase energy relation discussed 
in Section 2, these data are consistent with a large helicity-flip pomeron coupling. The real 
and imaginary parts of r 5 cannot be separately determined from the measurement of An at 
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this one value of t , but they could both be determined by measuring the t-dependence at 
RHIC because the deviation from the pure CNI shape is extremely sensitive to Kerr,. 



3.3 Iso-scalar part of the helicity-flip from ir ± p scattering 



Detailed Regge fits were made to spin dependent measurements in the 1970's |57L 



5q , p9j , p0| . At the low energies at which those measurements were made, there were quite 
large asymmetries observed. It was found that these were mainly due to the low-lying Regge 
trajectories and were not very sensitive to the pomeron couplings. The parameters that were 
found do predict a very small (< 10%) ratio of the flip to non-flip residues for the pomeron, 
but the parameters are uncertain because of this insensitivity of the fits. 

Polarization in up elastic scattering at high energy is mostly due to interference of the 
pomeron non-flip amplitude with the helicity-flip part of the p-Reggeon. As a consequence, 
the polarization has different signs and is nearly symmetric in ir^p scattering. It decreases 
with energy as 

^(M)oe(^J , (30) 

where a p (t) w 0.5 + 0.9 t and ap(t) w 1.1 + 0.25 t. The polarization has a double-zero 
behavior at t ~ —0.6 GeV 2 , which is correlated to the change of sign of a p (t) at this point; 
see Fig. 3. This effect can be understood as a result of destructive interference with the 
pCEDiP cut. An alternative explanation involves the wrong signature nonsense zero |58| (zeros 
in the residue and in the signature factor of the p-reggeon). 

At very high energies this part of the polarization vanishes, and one can hope to 
detect an energy-independent contribution of the pomeron. Unfortunately, available data 
are not sufficiently precise yet. One can eliminate the large background from the p <S> JP 
contribution by adding the data on polarization in ir^p elastic scattering, 

Z np (s,t) = 6 + (s,t) A* N +p (s,t) + 6-(s,t) A* N ~ p (s,t), (31) 

where 

6 ± ( 8 , t) = + 2<1 P(M) , (32) 
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and 



<£*>(., t) = M M_ exp ( B fp ( ) . (33 ) 



Therefore, Eq. (|32|) can be rewritten as, 

6+(s,t) 



2 7(M) 



^ TTtW' (34) 

where 

7 (, : 0-1 I <~^'. (30) 



a 



tot 



The difference between the elastic slopes AB = p — 5^ p is related to the position of 
the cross-over point to ~ —0.15 GeV 2 , which is nearly energy independent |61] in this energy 
range since 



a 



IT pi 

tot 



AB(s) = — In "SsTT ' ( 36 ) 



To find 5±(s, £) we fit the data on a^ ot p (s) |Q with the expression 

/ s \«p(o)-i /s\ a /(°H T / s \«p(o)-i 

We fixed a P (0) = 1.1, a/(0) = a p (0) = 0.5, s = 1 GeV 2 and found a P = 12.4 ± 
0.03 mb, o f = 40.8 ± 0.26 mb, a p = 5.1 ± 0.07 mb. 

Due to the cancellation of the isovector terms in Eq. fl3T|) E Tp (s, t) is dominated by 
the interference of the pomeron with the leading isoscalar reggeons. In up scattering this 
can only involve /-reggeon interference. 

Since the main part of the polarization cancels in the sum, the data have to have 
sufficiently high statistics in order to use Eq. (|3T|). This is why we could only use the 
low-energy data at momenta pi = 6 — 14 GeV/c, depicted in Fig. 3. 
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Figure 3: Polarization in up elastic scattering in the energy range 6 — 14 GeV . The data 
points are from ^/-/|6^/. The curves show our fit with the parameterization Eq. ( fEp - Eq. 



We performed a fit of A^ N p (s,t) with the parameterization 



A *±, ,s _ Z np (s,t)±A 7rp (s,t) 
Av(Mj ~ 2 5 ± (s,t) 



(38) 



where 



\t\ „ / g \ "flW-af (*) 

K P {s, t) = — e fll * (a 2 + ogt) (t - a 4 ) 2 K(s) - 



and 



Itl / S \ a R (t)-a P (t) 

S Tp (s, t) = - — e ait (a 5 + a 6 t) K(s) - 



(39) 



(40) 



Here t is in (GeV/c) 2 . We use the same energy dependence for E^s, t) and A 7rp (s,t) 
assuming that a nit) = a p (t) = a fit). The factor 



a r I S \ a R( t )~ a P (*) 

^-i( s) = 1 + 2^1 ( — 



0> VS 



+ 



sm 



fa/(0) Vso 



tts(')- a p(') 



(41) 
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takes into account the contribution of the /-reggeon to the differential cross section. The 
parameter a\ corresponds to the difference between the slopes of the pomeron and /-reggeon 
amplitudes. The factor (t — a^) 2 is introduced to reproduce the double- zero behavior of the 
polarization clearly seen in data |63|]-|65|]. It is usually related to presence of an additional 
zero in the p-reggeon residue which is dictated by duality at a p = 



The result of the fit is shown by the solid curves in Fig. 3, and the values of the 
parameters collected in Table 2. 





a 2 


a 3 


(24 


a 5 


a 6 


2.5 
±.2 


12.5 
±1.0 


68.5 
±8.2 


0.7 
0.01 


0.77 
±0.13 


2.33 
0.59 



Table 2: Fitted values of the parameters ai 
The ratio of the helicity-flip to non-flip isosinglet amplitudes can be extracted from 



r f (t) - r P (t) = \ f« (05 + aet)^ tan (^^) 



(42) 



Here rp^j\t\/m 2 and rfy\t\/m 2 denote the ratio of the helicity-flip to non-flip amplitude 

corresponding to IP and / exchange, respectively. We assume here that the helicity-flip and 
non-flip amplitudes have the same phase which corresponds to dominance of Regge poles. 
We neglect the real part of the pomeron amplitude. If we assume factorization, then for pp 
scattering, asymptotically, Imrs(s,t) «rp(t). 



Thus, the combination Eq. fl42|) of spin-flip to non-flip ratios for iso-singlet amplitudes, 
which is quite difficult to measure directly, are fixed by this analysis with a good accuracy. 



r/(0) -r F (0) = 0.06 ±0.01 . 



(43) 



Unfortunately, without further information regarding 77, this does not restrict rp. The 
approximation of /-dominance of the pomeron yields rp = Tf, which obviously contradicts 
Eq. (|43D . The pion exchange model in Section 5.2 predicts values for both rp and r/ which 
are in pretty good agreement with Eq. ((43]). The Regge fits of |58| use 77 = 0. This would 
give rp = —0.06, a very interesting value as we will see in Section 5. However, this should 
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probably be disregarded because the fits of |58[ also set the pomeron helicity-flip coupling to 



zero. The fits of [^] give, assuming exchange degeneracy and using the u> Regge residues, 
Tf = 0.95/10.6 ~ 0.09 and so, from Eq. (pE3"|), Tp = 0.03. Since this result requires some 



theory that is not tested to this precision, this can be taken as provisional but suggestive. 

Another source of information on the isoscalar exchanges is pA scattering. This 
requires special attention which we leave to another occasion. 



4 Model-independent bounds and the energy depen- 
dence of helicity-flip 



The magnitude of rs depends on the scale 1/m chosen in Eq. fl2~T|), where m denotes the 
nucleon mass. This scale has been used conventionally for many years; it was probably 
chosen in analogy to the form of the one-photon exchange helicity-flip amplitude. It is not 
at all certain that this is the appropriate scale for the scattering of strongly-interacting 
particles with structure. It might be more natural for the scale to be set by the slope of the 



diffraction peak; i.e. the effective radius of the proton R(s) = y2£>(s), (we take this to be 

the definition of the quantity R(s), see Eq. (Q) and Eq. (|5*E| ) below.) Since this is a good deal 
larger than 1/m, the "natural" size of r§ might be expected to be larger than 1. Furthermore, 
it might very well be expected to increase slowly with energy, corresponding to the growth 
in the effective radius of the proton. This, of course, flies in the face of conventional wisdom; 
see the discussion of Section 3. 

It is natural to investigate if there is a theoretical argument that r 5 — > as s — > oo. 
We begin by remarking that for the pure Coulomb amplitudes this is not true, so we ask 
if there is something different about the hadronic amplitudes. One obvious difference is 
that experimentally <fi + grows faster than s, and presumably will eventually grow as s In 2 s, 
the maximum rate allowed by the Froissart-Martin bound [ftOfl . Let us see what the same 
arguments used to derive that bound yield when applied to 05. The partial wave expansion 



22 



for 6?, is [EE] 



[s,t) = J2(2J+l)fj( S )d J w (6) 
j 



E(« + D^/JW^l'tcosS), (44) 

,(Z,m) , 



where t = — 2fc 2 (l — cos#) and Pj' (cos 6) denotes the Jacobi polynomial in cos# | 6"7fl . From 
this one finds that 



Us,0) = ^E(2^+ ^^mP^il), (45) 



where 05 = m/v— t 05. Partial wave unitarity requires that [68 

2|//( S )| 2 < Im#( S )(l - lmf^(s)) < 1/4. (46) 

If we assume that this bound is saturated out to some L max (s) ~ kR(s), where k « \fs/2 

is the CM momentum, then using Pjl\ (1) = J (to be compared with -Pj(l) = 1 for the 

Legendre polynomials) we find that for s — > oo, 5 (s,O) goes as msR 3 (s) while 0+ goes as 

sR 2 (s), and so the natural scale for 5 (s,O) zs P(s), not 1/m. This means that unitarity 
and other general principles allow r§ to grow with energy; if the Froissart bound is saturated 



I mH ~ ^/i In s and r 5 ~ In s is allowed. Note that if L max ~ Vslns then a tot will grow only 
as Ins as favored by Block et al [25|; in that case, r 5 ~ In 1 / 2 s is allowed. 

The above argument assumed the same L max for 0+ and 05. This can, in fact, be 

proved as follows: one can bound Pj 1 ' 1 '* from below, parallel to Martin's argument for p, 
the Legendre polynomial, in the unphysical region | cos^] > 1. One then applies the same 
reasoning as he used for + to 05. The representation 



D (i,i)/ s 2 ^ 



2 J f 71 

— / d6 (x + \fx 2 - 1 cos 0) J_1 sin 2 (47) 

71 JO 



allows one to show that PjZi(x) ~ x J ' j\fl as J —>■ oo for x > 1. This is the same as 
the asymptotic behavior obtained for the Pi(x) by Martin, and so polynomial boundedness 
implies the same L max for 5 and + . 
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Notice that the same arguments applied to the double-flip amplitudes 4>±{s, t) or 
d(p2(s, t)/dt\ t =o will yield a natural scale of R 2 (s) and, correspondingly, a possible growth 



One can easily see that </>5 can grow faster with s than <p + without violating unitarity 

because of the factor of \/~—t. It is, naturally, an interesting question to determine to what 
degree the helicity-flip amplitudes saturate unitarity, even at energies where the Froissart 
bound is not saturated. Techniques using unitarity and partial wave expansions have been 
used in the past at low energy to obtain bounds on the helicity-flip amplitude in terms of 
cr tot , cr cl and B |J^, [70, [71] ; these bounds are comparable in size to mR/2, i.e. r 5 is found to 
lie between 2 and 3. 

We can make the discussion of unitarity more quantitative by transforming the scat- 
tering amplitudes to the impact parameter representation. To keep the discussion as simple 
as we can, let us do this for scattering of a proton on a spin target; as a 2 x 2 matrix, the 
scattering amplitude has the form 



The two-dimensional Fourier transforms of these into impact parameter space yields 
the profile functions gi(b, s) and g2{b, s): 



with energy as fast as In 4 s. 





where q = k' — k , q = \q\ and q 2 



—t for elastic scattering. 




(49) 



where 





(50) 



With this normalization 




(51) 
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and unitarity imposes, for each value of b, the condition 

2Im<f 1 (M) > \g~i(b,s)\ 2 + \g~ 2 (b,s)\ 2 . 



(52) 



(This equation is, in general, only approximate in b space, but it can be derived from the 



analogous partial wave inequality |68| if only the elastic scattering amplitudes are sufficiently 
peaked in t.) The bounds discussed earlier correspond to a uniform distribution in b for both 
amplitudes for b < R = L max /k. If this ^-distribution is translated into the t-dependence 
of the amplitudes near t — it implies that the slope of g 2 j 'y/—i is less than the slope 
B = L 2 nax /2k 2 of gi; in fact it is 35/5. 

A more conventional assumption is that the slopes of g± and g 2 j \f—i are the same. 
If, in fact, g 2 (s, q) = A (q/m) gi(s, q), with A independent of t then 



g~2{b,s) 



A dg^b, s) 



rn 



db 



(53) 



This is true in the optical model or in any other model where the potential shape or matter 
distribution is the same for spin-orbit force as for the purely central force. Then g(s, b) will 
be more peripheral than for the bound just discussed. It has nothing intrinsically to do 
with unitarity or saturation of the Froissart bound, and it is clearly interesting to determine 
whether it is true or not. 



The unitarity condition Eq. fl52]) imposes a bound on |A|, and the closer g~i(b, s) is to 
saturating unitarity, the stronger this bound will be. Approximating the t-dependence of the 
amplitude by a logarithmically shrinking diffraction peak and neglecting its real part gives 



9i(b, s) 



ia\s) 



2irR 2 (s] 



exp 



R 2 (s) 



(54) 



and 



9~2{b,s) 



2 ibXa(s) 



exp 



R 2 (s) 



(55) 



27rmR i (s) 

where here and in the rest of this section the energy dependent Regge radius of interaction 
is 



R 2 (s) = Rl + ±a' ]P \Yi (—) 



(56) 
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With this form for the amplitudes a(s) = 0t o t(s) via Eq. ([)!]). This will change at the 
next stage of the calculation. Here one finds, numerically, over a wide range of values of 
a(s)/2nR 2 ( y s) < 1 that 



IAI < mR t 



2irR 2 (s) 



(57) 



is required in order to satisfy Eq. (^2|). 

If a(s) grows faster than R 2 (s) with s as s — > oo, say as s Ap [31], then the amplitude 
Eq. fl5"4|) will eventually violate the unitarity condition Eq. (|52|) and the form must be modi- 
fied. It is well-known that the total pp cross section at available energy is still far below the 
Froissart-Martin bound; however, the bound Eq. (^) is already saturated at small impact 
parameters, even ignoring the helicity-flip piece |34], |72|| . In principle, unitarity is restored 
after all the Regge cuts generated by multi-pomeron exchanges are added [|35[| . A standard 
way of unitarization of the non-flip part of the pole amplitude ]58| is eikonalization; however, 
the presence of the helicity-flip component may lead to problems with unitarity. Indeed, an 
even number of repeating helicity-flip amplitudes contribute to the non-flip part, but all of 
them grow as a power of energy and have the same sign. Therefore, eikonalization of the 
helicity-flip amplitude alone does not save unitarity, which can be restored only after the 
absorptive corrections due to initial/final state spin non-flip interactions are included. The 
resulting profile function reads, 



gr ik (b, s) = 1 - exp[i<7i(6, s)] + { 1 - cosh 



2i\b 
m R 2 (s] 



9i(b, s) 



exp[i g x (6, s)]. 



(5f 



The first two terms on the r.h.s. of this equation correspond to eikonalization of the non-flip 
part of Eq. (|48|) . They obey the unitarity bound at any s and b. In the extreme asymptotic 
region where a(s) in Eq. (Q) and Eq. (|55| ) is much greater than R 2 (s) then the 6-dependence 
has the form of a "black disk", i.e. gi eth (b, s) = 1 at b < R(s) and vanishes at b > R, where 



R 2 = Ajpln(-) R 2 (s) , 



(59) 
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if er(s) ~ (s/sq) Af . Likewise, 



g 2 eik = -sinh 



mR^ls) 



exp[z#i(fc,s)]. (60) 



Problems with unitarity at b < R(s) may arise from the last term in Eq. (|58|). The 
condition Eq. (|5^) is satisfied if 

. „ , . mR 2 (s) , . 

|ReA| < — 26 • (61) 

The minimal bound corresponds to a maximal b = R, and s — > oo, 



a 



i 

I \ 2 



| Re A| < m J . (62) 

For reasonable values of a' P and A P we conclude that | ReA| < 1.6. This is not a severe 
restriction, and is valid only in the extreme asymptotic limit, beyond the RHIC range; 
numerical calculations give a much larger bound, of order mR at RHIC energies. 

Note that A is renormalized by the eikonalization process; the result, X eik (s) can be 
calculated numerically from 

(S) - 2Jdbbg^(b,s)- (63) 
Likewise, the total cross section will be modified from the input values a(s) and is given by 

vtot(s) =4vr Jbdblmgr ik (b). (64) 

These last two equations will have to be used for comparison with data. 

5 Models for the pomeron helicity-flip 

An early attempt to understand the spin structure of the pomeron coupling was made by 



Landshoff and Polkinghorne [73|. This model preceded the formulation of QCD, but used 
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some of its features in a model they called the dual quark-parton model. They argued 
that the t-dependence of the pomeron coupling was determined by the electromagnetic form 
factors of the proton and neutron. This led to the conclusion that the helicity-flip coupling 
is given by the isoscalar anomalous magnetic moment of the nucleons; in our notation r 5 = 
— 1 + /i n )/2 = —0.06. This relation has subsequently been obtained or conjectured 
independently in a variety of models based on QCD. The result is, however, model-dependent 
as we will see. 

5.1 Perturbative QCD 

There is a widespread prejudice that the perturbative pomeron does not flip helicity. 
It is true that the perturbative pomeron couples to a hadron through two t-channel gluons, 
and that the quark-gluon vertex u^^Uq conserves helicity. However, one cannot jump to 
the conclusion that the same is true for a proton. In QED the fundamental vertex has the 
same form but radiative corrections induce helicity-flip via an anomalous magnetic moment. 



Ryskin [|74]] evaluated the pomeron helicity-flip coupling by analogy to the isoscalar anoma- 
lous magnetic moment of the nucleon. Applying this analogy to the quark gluon vertex he 
found the anomalous color magnetic moment of the quark. Thus the quark-gluon vertex 
does not conserve helicity and one can calculate the helicity-flip part of the pomeron-proton 
vertex. Using the two-gluon model for the pomeron and the nonrelativistic constituent quark 
model for the nucleon he found [ |73fl 

Imr 5 = 0.13, (65) 

independent of energy. In the above one needs to introduce an effective gluon mass and if one 
takes a large effective gluon mass, m g « 0.75 GeV, this estimate is substantially reduced. A 



need for a large gluon mass follows from lattice QCD calculations |75[ and the smallness of 
the triple-pomeron coupling |76| . 



The spin-flip part of the three-gluon odderon was also estimated in j74] and the 
helicity-flip component was found to be nearly the same as for the pomeron. If this is so 
then the odderon-pomeron interference contribution to An vanishes. See Eq. (P) and Table 
3 in Section 6. 

An alternative approach is to note that helicity is defined relative to the direction 
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of the proton momentum, while the quark momenta are oriented differently. Therefore, 
the proton helicity may be different from the sum of the quark helicities 0. The results 
of perturbative QCD calculations show that the helicity-flip amplitude in elastic proton 
scattering very much correlates with the quark wave function of the proton. Spin effects 
turn out to cancel out if the spatial distribution of the constituent quarks in the proton 
is symmetric || 177]]. However, if a quark configuration containing a compact diquark (ud) 
dominates the proton wave function, the pomeron helicity-flip part is nonzero. The more 
the proton wave function is asymmetric, i.e. the smaller the diquark is, the larger is Imrs 
[@j IZZI" ^ s va l ue i n t ne CNI region of transverse momentum ranges from —0.05 to —0.1 and 
even to —0.15 for the diquark diameters 0.5, 0.3 and 0.2 fm, respectively. The commonly 
accepted diquark size is 0.3 — 0.4 fm; therefore, we conclude that |Imr 5 | does not exceed 
10%. 

Note that there is a principal difference in sensitivity to the proton wave function 
between the helicity-flip and the non-flip components of the pomeron. The former probes the 
shortest interquark distances in the proton (diquark), but the latter is sensitive to the largest 
quark separation (due to color screening). Correspondingly, the virtuality of the gluons in the 
pomeron is higher in the helicity-flip component since these gluons must resolve the diquark 
structure. This fact may be considered as a justification for perturbative calculations for the 
helicity-flip part, while their validity for the non-flip part is questionable. 

High gluon virtuality in the helicity-flip pomeron leads to a steep energy dependence. 
A prominent experimental observation at HERA is that the steepness of growth with energy 
of the total virtual photoabsorption cross section correlates with the photon virtuality Q 2 , 
i.e. with the qq separation in the hadronic fluctuation of the photon. Analyses of the data for 
the proton structure function ^(x, Q 2 ) performed in [78 shows that for a quark separation of 
the order of the mean diquark diameter one should expect the energy dependence ~ (s/so) ' 2 . 
This should be compared with the well known energy dependence of the non-flip amplitude, 
~ (s/sq) oa . Therefore if the perturbative QCD model is meaningful in this region we expect 
a negative Imr 5 with energy dependence (s/so) 01 . This prediction can be tested in future 
polarization experiments at RHIC whose energy ranges from s ~ 50 GeV 2 (with a fixed 
target) up to 25 x 10 4 GeV 2 . Imr 5 is expected to double its value in this interval. 
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Eventually this growth will cause the bound Eq. ( |62| ) to be violated. This occurs 
only at very high energy, well above the LHC energy, and so it is not important for our 
considerations. Nevertheless, it would be interesting to develop an eikonalization method 
that would lead to consistent unitary amplitudes. We believe the eikonalization procedure 
developed in j?9| is the appropriate technique. When the elastic amplitude depends on 
transverse separation between partons, as it does here, the measured amplitude is the result 
of averaging over different transverse configurations: 

/(M) = </(M,V0>*, (66) 

where if) characterizes the transverse configuration and the averaging is weighted by the 
probability to be in configuration if). Correspondingly, eikonalization has to be done first for 
a given configuration if) and only then averaged : 

r k (b,s) = (r k (b,s,ip))^. (67) 

For a given partonic configuration if) the energy dependence of the helicity-flip and non-flip 
components must be the same since, as stated above, it depends only on the transverse 
separations. Therefore, restriction Eq. (^) applies except that the pomeron intercept de- 
pends on if} and unitarity is satisfied for each if). However, the weight factors are different 
for the helicity-flip and non-flip amplitudes and the averaging results in a higher effective 
intercept for the helicity-flip component. The detailed predictions of this procedure remain 
to be worked out. 

5.2 Pion exchange model 

A nucleon is known to have a pion cloud of large radius. Since the helicity-flip am- 
plitude is proportional to impact parameter, it is natural that a substantial fraction comes 
from inelastic interaction of the projectile hadron with virtual peripheral pions. This contri- 
bution is related through the unitarity relation to a pomeron-nucleon vertex (in the elastic 
hadron-nucleon amplitude) shown in Fig. |j. It is known that the main contribution to the 
pion cloud comes from the virtual transitions iV — > ttN and iV — > 7rA, which corresponds to 
the two graphs depicted in Fig. |j. This model for the pomeron-nucleon coupling was sug- 
gested in PD[ . They predicted Imr 5 « 0.016 (lns) 3//2 . This quite a steep energy dependence 
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Figure 4: Pomeron coupling to a nucleoli via two pion exchange. 

originates from the radius of the pion cloud which is assumed to be proportional to y/hi s. 
A more detailed analyses was undertaken in |yj. An interesting observation of this paper is 
a strong correlation of the value of r 5 with isospin in the t-channel. Namely, for an isoscalar 
exchange (IP, /-reggeon) the two graphs in Fig. [| essentially cancel in the helicity-flip, but 
they add up in the non-flip amplitude. It is vice versa for an isovector exchange (p-reggeon). 
This conclusion is consistent with Regge phenomenological analyses of experimental data 
(see e.g. f82j). 

In order to fix the parameters of the model a detailed analysis of data on inclusive 
nucleon (pp — > p(n)X) and A (pp — ► A ++ X and 7i + p — ► A ++ X) production was performed 
in [plj . These reactions correspond to the unitarity cut of the graphs in Fig. [1. The calcula- 



tions in EI] led to a positive value of Imr 5 = 0.06 for the pomeron (0.15 for the /-reggeon). 



This nonperturbative contribution has the opposite sign to what follows from perturbative 
calculations and may partially compensate it (see discussion in ||). 



5.3 Impact picture 

An impact picture approach, which was derived several years ago [|83], 0, |85| , describes 
successfully pp and pp elastic scattering up to ISR energies . It led to predictions at very 
high energy, so far in excellent agreement with the data from the CERN SPS collider and 
the FNAL Tevatron and others, which remain to be checked at the Large Hadron Collider 
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under construction at CERN . The spin-independent amplitude reads at high energies 



impact,* = ls i j Q ( h J-ft(l _ e -lW))bdb , (68) 
JO 

where the opaqueness Qq, which is assumed to factorize as £Iq(s, b) = So(s)F(b 2 ), is associated 
with the pomeron exchange. The energy dependence is given by the crossing symmetric 
function 

S {s) = s c / ln c ' s + u c j ln c ' u, (69) 

which comes from the high energy behavior of quantum field theory. In Sq(s) above, u is 
the third Mandelstam variable and both s and u are expressed in GeV 2 . Note that So(s) 
is complex because u is negative. The phenomenological analysis leads to the values of the 
two free parameters c = 0.167, c' = 0.748 and the real part of <ft l ™ pact (s, t) results from the 

phase of S (s). The t-dependence of 0![ npoc *(s, t) is driven by F(b 2 ), which is related to the 
Fourier transform of the electromagnetic proton form factor and, as a result of a simple 
parametrization which can be found in | 83fl , F(b 2 ) is fully determined in terms of only four 



additional parameters. 

The spin structure of the model was also studied and it allows a rather good descrip- 
tion of the polarization data, up to the highest available energy, i.e. pi = 300 GeV/c |86| . 
At the RHIC energies, the spin dependent amplitude reads 



^ ac \s,t) = is / J 1 (bV^i)n 1 (s,b)e- n ° {s '%db , (70) 



o 

where b) is the spin dependent opaqueness, corresponding to the helicity-flip component 
of the pomeron. It also factorizes as Q\(s,b) = Si(s)F s (b 2 ), where Si(s) is obtained from 
S (s) and we have 

Si(s) = -5-(c-cyiiis) + (s->i0 . (71) 
In s 

F s (b 2 ) is simply related to F(b 2 ) according to F s (b 2 ) = bu(b 2 )F(b 2 ), where u(b 2 ) is a smooth 
function which is not very precisely known. The important point is its value ujq for very 
small b and by fitting the data, it was found that ujq = 0.06 GeV. This leads to a value 
Imr5 ~ —0.06, if one assumes that the flip component of the pomeron is normalized at 
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t = 0, by the nucleoli isoscalar magnetic moment ||TT|. This is at variance with the exact 
results one obtains in the impact picture, which are shown in Fig. [5] at two different energies. 
It is interesting to remark that Imrs(t) increases with energy, in a way pretty much consistent 
with what was mentioned above in Section 4. 
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Figure 5: Imrs(t) calculated in the impact picture for two energy values: = 50 GeV 
(dashed curve) and y/s = 500 GeV (solid curve) . 



6 P-independent determination of 02, </>_, (f>4 and 05 

In this section we would like to demonstrate that, in principle, by making use of both CNI and 
hadronic interference at small t it is possible to determine all the spin dependent amplitudes 
at t = independent of knowledge of the beam polarization P\ and P2 provided only that 
they are stable and non-zero. This is very interesting, perhaps surprising, in its own right. If 
it proves to be practical, it would permit the use of elastic pp scattering as a self-calibrating 
polarimeter. It is important to emphasize right at the beginning — we will not repeat this 
every time the issue occurs — that the method involves several ratios of very small quantities; 
the precision required to do this may be beyond the reach of practical experiment at this 
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time. However, very little is known about the amplitudes now, and so we cannot evaluate 
this. Even if the complete process we describe cannot be carried through, much of what 
follows should be useful in constraining the amplitudes at small t. 

The method requires the use of asymmetries with both longitudinal and transversely 
polarized beams; it will not work unless data are available with both configurations. Here 
we work only to order a and so only amplitudes that are large compared to the next or- 
der correction can be determined from the formulas given below. This could probably be 
improved upon if necessary; at the present, experiment will probably not be able to probe 
amplitudes below that size and so we have not pressed on in this direction. We assume that 
the polarized beams have the same degree of polarization P in either configuration; since 
they are produced from the same initial configuration by rotation this is almost certainly 
true. For simplicity in writing we assume both beams to have the same polarization; this 
may very well not be true but it is trivial to correct the formulas for this. 

We work with the experimentally measured asymmetries which are given by PA N , 
P 2 Ann, etc. These will contain singular terms as t — > coming from the interference 
between the one-photon exchange and the hadronic amplitudes. To order a the asymmetries 
Ann, Ass and All are singular as 1/t and An,Asl are singular as \jyf—i. So we write 



m^J-i da a tot 

PA N — = -aa N + — b N t + 

a tot at 87T 

* o2 a da i a tot 



P All -77 = aa L L + -77— b L L t + 



a tot dt 8n 

1 r-,2 a da t a tot 



P Ann —r — cnaNN~\~~z — bNNt~\~. 



0tot dt 8n 

t 2 a da , fTtot 



P A ss -r = aa ss + —b S st + ■ ■ ■ 



a tot dt 87r 

my^i 2 A da i cr tot 



PAsl-77 = aasL + ^bs L t + .... (72) 

In the following table we give expressions for the various which we sometimes refer to 
as the enhanced pieces, and bi which we refer to as the hadronic piece. The a^s are linear 
in the hadronic amplitudes while the bi are bilinear. Here we omit terms of order at which 
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are small; we will return to show how this can be corrected for, if necessary. Note that the 
usual exponential t-dependence of the hadronic amplitudes will enter only at order t 2 or at; 
likewise, for the quantities <\)\ and 04. In this approximation Ann = Ass- For notation see 
Section 2. 



asymmetry 


a. 


h 


A N 


P{f(l + / 2 )-/ 5 } 


PI 5 {(p- p 5 ) + / 2 (p 2 - P 5 )} 


Ann 


P 2 R 2 


p 2 R2{ P + j- 2 +i-(p- + j- 2 )} 


A ll 




pi{R_(-L + p ) + R 2 I 2 (-L + P2 )} 


A SL 


P 2 {I(R^ + R 2 )} 


P%{(I_ + I 2 )+p 5 (R_ + R 2 )} 



Table 3: The first two terms in an expansion in t of the various asymmetries. 

(The possibility of using the electromagnetic and hadronic pieces of A N n and A L l to 
determine the real and imaginary parts of 2 and 0_, when the polarization is independently 
known was noticed in 

We also need the cross section differences 

Acr T = -2/ 2 o- tot , 

Aa L = 2I_a tot . (73) 

Fits to the data will determine ciN,bN,a>NN,bNN, etc. The strategy will be to take ratios 
of two quantities that are either linear or bilinear in the polarizations to obtain ratios of 
amplitudes which will then be independent of the polarization. We will find that there are 
enough of these ratios to solve for all the amplitudes, provided that at least one of 2 and 
0_ is non-zero. Indeed, the system is over-constrained, and the procedure we describe here 
is not unique. We carry it through here to demonstrate that a solution exists; the optimal 
method will no doubt depend on the experimental situation. If both 2 and 0_ turn out to 
be unmeasurably small, the method fails at step one. 

Let us begin with the ratios of the four measured asymmetries: the total cross section 
differences and the enhanced parts of Ann and All- From these one can get immediately 
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the ratios of real to imaginary parts for 02 and 0_ 



p 2 = -2a NN a tot /P 2 Aa T , 

p_ = 2a LL a tot /P 2 Aa L . (74) 

This fixes the phase of both amplitudes 02 and 0_. From the same four measurements a 
third independent ratio can be formed; either 

P 2 Ao"l ~~ I { } 



(76) 



or 

O-LL R- 
dNN R% 

will fix the ratios of the magnitudes of 02 and 0_. We will use the latter in the following. 

In order to completely fix the magnitudes, one more ratio is needed. Either b^N/dNN 
or bLi/cLLL will do. Examination of the table will reveal that either of these quantities 
depends only on I2 or, equivalently, i_ in addition to the ratios just determined; the unknown 
I2, say, is thereby related linearly to the ratio with known coefficients: 

r _ aiL ((Pll/oll) - l/p--p) 

1 2 — — 7-1 • {<<) 

dNN 92 + I/P2 

At this point, one has enough information to determine the polarization because one 
can calculate R 2 , 1- and R_ from Eq. ( [77D and the previously determined quantities: one 
uses either o,ll or cl^n in 

P 2 = ?™ (78) 
R2 



or 

t-,2 a LL 



R. 

to obtain 



(79) 



a 2 NN + {P 2 Aa T /2a tot y 



b LL -P*(Aa L /2a tot )-pa LL 1 ] 
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This equation is valid in all the various degenerate limits except the case 02 = 0, both real 
and imaginary parts, in which case it is indeterminate and one must work harder. 

Barring this exceptional case, one is in principle done because with this P — presumably 
the sign ambiguity will not present a problem — one can use the table to calculate 05 from 
An] Asl is not needed. To give an idea of the sensitivity of Ann to R2, the curve for Ann 
has essentially the same shape as the CNI curve for An and, for R2 = 0.02, the height at 
the maximum is about 2%. It may very well happen that Ann is measurable but that the 
error is too large for this to provide a precision measurement; ±1% would not be useful in 
the example just cited. Here, too, one may benefit from pressing on: an error of ±1% in I 5 
would be far better than is required because it is applied to a term of order 1 in A N . 

Going further requires bringing in 6/v / cln and bsi/asL- Each of these can be used 
to express I 5 in terms of p 5 and measured quantities. By equating these two expressions an 
equation for p 5 is obtained. The result is 

1 / (W«SL)((W«jv) +P + R2) {b N /a N ){I 2 + I-) \ 

P5 b N /aN + bsL/asL\ l + h R2 + R- J ' 1 ' 

If this is then inserted into the equation for, say, b N /a N then J 5 is determined since we have 

J 5 K(b N /a N + b SL /a SL ) \ l + h R 2 + R-j' 

Notice that there are no quadratic ambiguities in any of these determinations. This is valid 
in all degenerate cases as well, as can be easly checked; it only fails if both 2 and 0_ vanish. 
The problem then becomes identical to that of a proton scattering off a spin particle for 
which one cannot calculate the spin dependence without knowing P. 

This procedure can be extended to apply to the case where the two spin 1/2 particles 
are distinguishable, as in p - 3 He scattering. The part concerning Ann, All and Ass is 
identical. There are two new quantities to determine, p§ and h, and there are two additional 
equations, effectively from bsi and clsl and from b' N and a' N . These can be solved just as in 



Eq. (|1D and Eq. (||). 



One can imagine a number of special cases. An interesting case is pure pomeron pole 
dominance. In that case (cf. Section 2) + , 02 and 05 are all in phase while 0_ = 0. In 
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this very simple case, which should be easily checked experimentally, bsi/bNN = h and so 
05 is determined in terms of measured quantities. Equivalently, one can use the ratio of the 
hadronic piece to the enhanced piece of A$l- The corresponding ratio of the hadronic piece 
of All to the enhanced piece of Ann determines I2 so everything is fixed: 

K b SL 1 



2a SL (p + l/py 

b LL 1 
a N N (p+l/p)' 



a) 



One can easily take into account the Bethe phase corrections to this procedure. Evi- 
dently, it will modify only the Oj and has no effect on the 6j. We have already seen in Section 
3 that because S is so small and because it enters An only by multiplying small quantities p 
and r 5 or r 2 , it can be safely neglected in to the accuracy that we are working. The correc- 
tions to Ann and All are very similar; so, cinn — ► P 2 (R2 + 8 h) and o,ll —> P 2 {R~ + 5 I_) 
to lowest order in 5; thus Eq. ([74]) becomes 

P2 + 5 = -2a N N cr t ot/-P 2 Acr T 

P- + 5 = 2a LL a tot /P 2 A(TL. (84) 

Since 5 is a known quantity, the values of P2 and p_ can be determined for use in the 
subsequent steps. 

We now return to the at corrections; these are small but they may need to be taken 
into account in order to use this method if the amplitudes <p-,<f)2 and <p 5 are quite small. 
The explicit expressions for these terms are given in detail for all of these asymmetries in || . 
There are several sources of these corrections. The most important arises from the slopes 
of the forward hadronic amplitudes, call them Bi. In the purely hadronic part they appear 
only in order t 2 but, through interference with the Coulomb singularities in either + or (j) 5 , 
they contribute to 6j. It is very likely that the slopes for the helicity-flip amplitudes are 
not very different from the non-flip B + , a factor of 2 at most; cf. the discussion in Section 
4. To the degree that they are the same the correction to 6; is just a47iB + /a to t- This 
corrects the corresponding ratio fri/dj by a known amount and can be simply accounted for. 
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To the degree that the slopes Bi are different this procedure leaves behind a correction of 
(Bi — B + )/2 multiplied by one of the presumably small amplitudes 02, 0- or 05 an d so is at 
a level of about 1CT 3 . The forward slopes of the Coulomb amplitudes can be taken account 
of, in exactly the same way. 

There is a correction to the real part of 2 equal to 2are 2 /4m 2 which is about 0.01; this 
can simply be added into R 2 and everything goes through as before. The term proportional 
to at arising from |0 5 | 2 is of order 10~ 3 and so can be ignored. 

Finally there is the heretofore unmentioned 04 which vanishes linearly with t as t — > 0. 
Although the amplitude never enters the enhanced piece, it does enter through interference 
into the linear term in t. One guesses that its contribution will be negligible, but since 
nothing is known about it, one would like to make sure that it can be controlled. Indeed, 
it can in principle be determined by this method: this amplitude can be removed from the 
first steps of the game by using (Ann + Ass)/2, instead of Ann- The determination of the 
amplitudes 0_ and 2 goes through as before. Then by considering (A N n — A ss )/1 one can 
determine R^. This can be used to correct A$l which can in turn be used to fix p 5 . Finally, 
then 6jv can be used to fix 7 4 and everything is determined. 

We don't want to oversell this method for self-calibrating CNI polarimetry; we realize 
it is experimentally very uncertain. However, even if the essential asymmetries are too small 
for this method to succeed, this linear parametrization, making use of the CNI enhance- 
ments, should prove useful for determining the amplitudes at t — 0, when the polarization 
is independently measured. Furthermore, we find it interesting that it is possible, at least 
in principle, to determine all of the spin dependent amplitudes without knowing the beam 
polarization independently. 

7 Conclusions 

Motivated by the need to have an accurate knowledge of the proton-proton single helicity- 
flip amplitude 05 at high energies, in order to devise an absolute polarimeter for use in the 
forthcoming RHIC spin program, we have examined the evidence for the existence of an 
asymptotic part of 5 which is not negligible compared to the largely imaginary average 
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non-flip amplitude + = \{(f>i + 03 ) at high energies. There is a general prejudice that 

r 5 = mfe/yf—i Im0 + will be negligibly small at high energies, say for p L > 200 GeV/c, and 
we have tried, using various techniques, to assess the validity of this belief. We have explained 
how certain characteristics of the dynamical mechanisms are linked to the behavior of the 
helicity amplitudes at high energies and small momentum transfers, namely their growth with 
energy, their phases, their small-t behaviour, and relations amongst them. On the basis of 
rigorous analytical methods we have demonstrated that the same fundamental assumptions 
which lead to the Froissart bound, |0 + | < s In 2 s , permit r 5 to grow like Ins. This surprising 
result implies that there is nothing in principle to stop 05 from remaining large, or even 
growing, relative to 0+ at high energies. However, other methods of analysis, based either 
on information at low to medium energies, or based upon dynamical models, do suggest a 
small 05 at RHIC energies, typically |r 5 | < 15%. 

Experimentally, for the region of interest to us, the best constraint on 5 comes from 
the measurement of A N in the CNI region at p L = 200 GeV/c. Assuming that the phase 
of 05 is the same as that of + — a sensible assumption for an asymptotically surviving 
contribution — one finds |r 5 | = 0.00 ±0.16. However, freeing the phase yields |r 5 | = 0.2 ±0.3 
and a phase difference between 5 and + of 0.15 ± 0.27 radians. We believe that the former 
value is the more reliable. There are conflicting non-perturbative estimates of r 5 at t = 0. By 
attempting to link helicity-flip to the isoscalar anomalous magnetic moment of the nucleon, 
Landshoff and Polkinghorne arrive at Imr 5 = —0.06 . This result is supported by an eikonal 
analysis of BourrelySoffer and collaborators, who find Imrs = —0.06 when the nucleon 
matter density is taken equal to the charge density. However, a more realistic choice of 
matter density leads to Imr 5 = -0.018 at y/s = 50 GeV and -0.026 at y/s = 500 GeV. 
Surprisingly, a study by the ITEP group, based upon the importance for helicity-flip of the 
peripheral interaction with the pion cloud in the nucleon, and which should therefore not 
be too different from analyses based upon the matter density, yields Imr 5 = 0.06 i.e., of 
opposite sign to the above mentioned results. On the other hand Ryskin has attempted 
to calculate the anomalous colour magnetic moment of a quark, based upon a mixture of 
perturbative QCD and the constituent quark model, and linking the result to 05 obtains 
Imr 5 = 0.13 i.e., of opposite sign to the results based upon the electromagnetic anomalous 
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moment. Perturbative QCD attempts by Kopeliovich and Zakharov to link the existence of 
helicity-flip to the transverse momentum of the consituents turn out to be very sensitive to 
the form of the nucleon wave-function. If the wave function contains a significant component 
corresponding to a compact scalar (ud) diquark they find that Im r 5 increases in magnitude as 
the diquark size D decreases. Quantitatively Imr 5 = —0.05 — > —0.15 for D = 0.5 — > 0.2 fm. 

In summary while the various approaches give results which differ in sign and magni- 
tude, and while it is not clear to what extent perturbative and non-perturbative approaches 
overlap, it seems reasonable to assert that |r 5 | < 10% at RHIC enrgies. This level of ac- 
curacy is unfortunately inadequate for the needs of an absolute polarimeter. We have also 
studied the amplitudes 0_ = \{4>i — 4>^) and 2 . There is persuasive evidence both from 
experiment and from dynamical arguments that 0_ is exceedingly small at high energies: 
\<j)-/<j)+\ < 10~ 3 for energies beyond pl = 200 GeV/c. The case of 2 is less clearcut. There 
is experimental evidence, but from relatively low energy measurements of Aa T , that Im0 2 
drops from —6% — > —0.4% for p^ = 2 — > 6 GeV/c. And there is evidence from charge 
exchange scattering that the 1=1 part of 2 is very small at higher energies: |r 2 | < 0.006 
at pl = 270 GeV/c. On dynamical grounds we expect |r 2 | — * 0, but the argument is not 
conclusive. 

Finally, we have demonstrated the surprising result that proton-proton elastic scat- 
tering is self analysing, in the sense that all the helicity amplitudes can be determined 
experimentally at very small momentum transfer, without a knowledge of the magnitude of 
the beam and target polarization. The experimental procedure for doing this is complex, 
but once carried out successfully it would permit the calibration of a CNI polarimeter which 
could then be used very simply for routine measurement of the beam polarization. 
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